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Abstract. I recall the development of quarkonium quantum mechanics after 
the discovery of T. I emphasize the empirical approach to determining the force 
between quarks from the properties of cc and bb bound states. I review the 
application of scaling laws, semiclassical methods, theorems and near-theorems, 
and inverse-scattering techniques. I look forward to the next quarkonium spec- 
troscopy in the Be system. 



PROLOGUE 

I am very happy to share in this celebration of the twentieth anniversary of 
the discovery of the 6-quark. The upsilon years were a very special time for the 
development of particle physics. Reviewing the events of two decades ago, I 
was struck not only by the pace of discovery, but by how easy it was, and how 
much fun we had. As a young physicist of that time, I am grateful not only 
for the excellent science, but also for the excellent people that quarkonium 
quantum mechanics gave me an opportunity to know, work with, and even 
compete with. They taught me a great deal about physics and life. 



THE EMPIRICAL APPROACH 

Charmonium quantum mechanics was already well-launched when the T 
came along. Appelquist & Politzer [1] had shown that nonrelativistic quantum 
mechanics should apply to QQ systems. The Cornell group [2] had shown 
the predictive power of the nonrelativistic potential-model approach using a 
"culturally determined" potential, 

^) Internet address: quigg@fnal.gov. 

^) Fermilab is operated by Universities Research Association Inc. under Contract No. DE- 
AC02-76CH03000 with the United States Department of Energy. 



V{r) = -- + -, 



(1) 



Eichten & Gottfried [3] had anticipated the spectroscopy of bb, predicting 

M{V) - M(T) « 420 MeV (2) 



|[M(^')-M(^)] 



All this meant that the deductive approach — assuming a form for the in- 
terquark potential and calculating the consequences — was in very good hands. 

Jon Rosner and I were motivated to take a complementary empirical ap- 
proach by the way the quarkonium problem came into our consciousness. The 
facts we had at our disposal in the summer of 1977 were these: 



E288 


M(T') - M(T) 


M(T") - M(T) 


Two-level fit 


650 ± 30 MeV 




Three-level fit 


610 ±40 MeV 


1000 ± 120 MeV 


M{ijj') - M{ip) 


590 MeV 





We were much impressed with the fact that the T' — T spacing is nearly 
the same as — ip. In an off-hand conversation with Bernie Margolis, who 
was also visiting Fermilab, Jon asked if he knew what kind of potential gave 
a level spacing independent of the mass of the bound particles. When Bernie 
said that it was probably some kind of power-law, Jon set off on the path 
that led to a potential V{r) = as e ^ 0, and eventually to the logarithmic 
potential. Accurate numerical calculations and a scaling argument showed us 
that the logarithmic potential, V{r) — Clog(r/ro) indeed gave a level spacing 
independent of mass [4] . 

The logarithmic potential gives a good account of the (then known) ip and 
T spectra, as shown in Figure 1, but is not unique in doing so. It was an easy 
matter to produce a modified Coulomb + linear potential that gave equal 
spacing for the ip and T families, but for no other quark masses. Now, the 
logarithmic potential is the solution to an idealized statement of the experi- 
mental facts. We expect it to be a good representation of the interaction in 
the region of space that governs the properties of the narrow ip and T states, 
but we have no reason to attach fundamental significance to it. It is mildly 
intriguing that a logarithmic confining interaction emerges from the light-front 
QCD approach when the Hamiltonian is computed to second order [5]. 



SCALING THE SCHRODINGER EQUATION 

The Schrodinger equation for the reduced radial wavef unction in a potential 

V{r) = Xr" is 



2pL 



E-\r 



2/xr2 



u{r) — . 



(3) 
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FIGURE 1. Level schemes of the ip, and hypothetical ^ families in (a) nature, (b) the 
logarithmic potential, and (c) a Coulomb + linear potential. (Prom Quigg & Rosner [4].) 



With the substitutions 



l/(2+u) 



E 



V2/i|A|y 



-2/(2+1^) 



(4) 
(5) 



and the identification w{p) = u{r), we can bring the Schrodinger equation to 
dimensionless form [6], 



w"{p) + 



e - sgn(A)p'' 



£{£+!) 



w{p) = . 



(6) 



The substitution (4) means that quantities with dimension of length scale as 
[L] oc whereas (5) tells us that quantities with dimensions of 

energy scale as [Ai?] cx (/i)~''/*^^"^'^)(|A|)^/(^+^). The scaling behavior in several 
familiar potentials is shown in Table 1. 



TABLE 1. How length and energy observables 
scale with coupling strength |A| and reduced 
mass fj, in power-law and logarithmic potentials. 



Potential [L] [E] 

Coulomb UAW^ JAW 

Log: y(r) = Clog r (C^ij-i/^ q^o 

Linear M"'/^|Ap/' 

SHO (m|a|)-i/* m"'/'|a|^/2 

Square well (m|A|)'^ 



Quarkonium Decays 

The scaling laws have immediate applications to the matrix elements that 
govern quarkonium decays. Electric and magnetic multipole matrix elements 
have dimensions 

{n'\Ej\n)r^L^, (7) 
(n'|Mj|n)~LV/x. (8) 

Radiative widths are given by 

r(Ej or Mj) ~ p'^^^+^\{n'\Ej or Mj|n)|^ , (9) 

so that transition rates scale with mass as 

r(Ei) (10) 

T{M1) ^ l^-^^'+^'^/^^+'l (11) 

Probability densities have dimensions L~^. Accordingly, the wave function 
squared at the origin scales as 

|*(0)|2 ~ , (12) 

so the leptonic width of a vector meson scales as 

r(V° ^ £+£-) = 167ra2eJ|^(0)|7M(V°)2 (13) 

Combining (10) and (11) with (13), we see that 

r(Ei)/r(£+r) ~ 12-^^+"^^^+"^ (i4) 

T{M1)/T{tr) ~ ^-3(i+-)/(2+<^). (15) 

For a potential less singular than a Coulomb potential, radiative decays be- 
come relatively less important than leptonic decays, as /j, increases. 
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FIGURE 2. Expectations for the leptonic widths of T and T'. The lower bounds (17) are 
indicated for et — —■^ (sohd lines) and Cb — ^ (dashed lines). The shaded region shows the 
widths predicted for e;, = — on the basis of twenty potentials from [20] that reproduce the 
ip and ip' positions and leptonic widths. The data point represents the 1978 DORIS results. 



Measuring the 6-quark's charge 

For a power-law potential with z/ < 1, the scaling law (12) implies that 

i*6(o)r>— i^e(o)r, (16) 

rric 

which leads to 

The inequality holds for more general potentials than power laws. We can 
prove it for the ground state for any monotonically increasing potential that 
is concave downward [7]. For excited states we have given a derivation for the 
same class of potentials in semiclassical approximation. 

The numerical bounds that follow from (17) are indicated in Figure 2. Mea- 
surements presented at the 1978 Tokyo Conference by the collaborations work- 
ing at the DORIS storage ring [8], 

T{T^i+i-) = 1.26 ± 0.21 keV , 
r(T'^£+r) = 0.36 ± 0.09 keV , 

ruled decisively in favor of the = — | assignment. The fifth quark was 
indeed bottom. 
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FIGURE 3. Comparison of "^Si levels of the ip and T families, circa 1979. 

The Order of Levels 

The discovery in 1975 of the narrow resonances Pc and x confirmed the 
theoretical expectation that the 2S and 2P levels are not degenerate in char- 
monium, as they would be in a pure Coulomb potential. In response to the 
question of what the 2S-2P splitting says about the interquark potential, 
Andre Martin and collaborators (R. Bertlmann, H. Grosse, J.-M. Richard, 
and others) constructed a series of elegant theorems on the order of levels in 
potentials [9,10]. 

After early data on the Upsilons confirmed the similarity of the J/t/j and T 
spectra, as shown in Figure 3, Martin used the scaling laws to deduce a simple 
power-law potential. 



V{r) = -8.064 GeV+ (6.898 GeV)(r • 1 GeV) 



0.1 



(18) 



The Martin potential [11] has served as a very useful template for the J/ip, T, 
and even (j){ss) families. In addition, it has led to many informative predictions 
for the masses of baryons containing charm and beauty. 



A Priority Dispute with Isaac Newton 

About three years ago, Jon Rosner telephoned me to say that Professor 
Chandrasekhar had just advised him that we were in a priority dispute with 
Isaac Newton. "Capitulate at once!" I said. "Newton can be very terrible." 

In writing his superb commentary on the Principia Mathematica [12], Chan- 
dra had found that Newton was the first to explore pairs of dual power-law 



potentials, and that he had mapped the Kcplcr problem into the harmonic 
oscillator. We had shown [13] — three centuries later — that the bound-state 
spectrum of an infinitely rising power-law potential, 



V{r) = Xr" (z/ > 0) , (19) 
is connected with that of a singular potential, 

V(r) = Ar^ - 2 < (P < 0) . (20) 
The paired Schrodinger equations for the two cases can be written as 



2fjL 
2pL 



u"(r) + 



v"{z) + 



E-Xz''- 



2/ir2 
2/iz2 



u{r) = , 
v{z) = , 



(21) 
(22) 



where {i^ + 2)(P + 2) = 4 and ^ = X{u/u)\ X = -E{u/uf, {£ + 1/2)^^ = 
(£+ 1/2^9^, andz = ri+^/2. 

The familiar quantum-mechanical correspondence between the Coulomb 
and harmonic oscillator problems emerges as a special case. For circular or- 
bits, Newton cites a relation between u and P equivalent to ours. We did 
capitulate, and so far we have not suffered any indignities at Newton's hands. 
And because Jon is a scholar, he and Aaron Grant have written an excellent 
historical review of the classical and quantum- mechanical analyses [14] . 



COUNTING NARROW LEVELS OF QUARKONIUM 

Eichten & Gottfried had argued that a QQ system with mg ^ rric would 
have at least three narrow ^Si levels [3]. As Kurt recalled in his talk, this 
implies a very rich spectroscopy, which figured prominently in the scientific 
case for CESR. However, the observed T' — T spacing is much larger than the 
420 MeV they predicted. Figure 4 shows that in the logarithmic potential, we 
would predict three or four narrow ^Si levels of T, in agreement with Eichten 
& Gottfried's expectation. This circumstance led us to ask how general is the 
expectation, and on what does it depend? 

Scmiclassical methods (whose power Taiji Yamanouchi had impressed on 
us) led us to a remarkable general result [15]: The number of narrow ^Si 
levels is 

The derivation is short enough to reproduce in full. 
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FIGURE 4. Plot of the Qq + Qq threshold relative to the (ground-state) l^Si QQ level 
as a function of the ratio of the heavy-quark mass mg to the charmed quark mass TOc for 
a logarithmic potential. Upper curve: reduced-mass and hyperfine corrections included. 
Lower curve (straight line): reduced-mass and hyperfine corrections ignored. Horizontal 
lines denote the 3'^Si, . . . QQ levels in this potential. (From Quigg and Rosner [15].) 



Set the zero of energy at 2mQ. The threshold for the dissociation of quarko- 
nium {QQ) Qq + Qq hes at an excitation energy 5{mQ). If V{r-) binds QQ 
states rising at least 6{mQ) above 2mQ, then the WKB quantization condition 
is 

/ dr[mQ{5{mQ)-V{r))Y/^^{n-\)'K , (24) 
Jo 

where tq is the point at which V^(ro) = 5 {mo). As Eichten & Gottfried had 
observed, Simq) = 2m(lowest Qq state) — 2mQ approaches a finite limit 6oo, 
independent of mq, as mq oo. This means that the only dependence of 
(24) on the heavy-quark mass is the explicit factor of ^Jvtiq on the left-hand 
side. We have, by inspection, the general result 

(n - i) oc ^ . (25) 

This universal behavior is realized in different ways for different potentials, 
as illustrated in Figure 5. The examples chosen are Vir) = r, V{r) = Inr, 
and V{r) = — r~^/^, for which AE oc fjP, and /i^^^^, respectively. All 

the levels fall deeper into the potential as the reduced mass fi increases, in 
conformity with the Feynman-Hellmann theorem. For the singular potential 
with z/ = — i, the levels spread apart as they sink into the well. For the linear 




potential, the levels are packed more densely as jj, rises. The logarithmic 
potential is an intermediate case in which the level spacing is independent of 
the mass and all levels fall into the well at a common rate given by Ei^jj!) — 
-I In (//'///). 

SEMICLASSICAL METHODS AND RESULTS 

The power of the WKB approximation for the counting problem encouraged 
us to explore other applications of semiclassical methods. Evaluating the 
nonrelativistic connection 

|.„(0)P^A(-)^ (20) 

in the semiclassical approximation, we connect the square of the s-wave wave 
function at the origin to the level density: 



l*n(0)| 



2 ^ M!^£;l/2^ 



47r2 



dn 



(27) 



(for a nonsingular potential) [16]. For example, in a linear potential V{r) = Xr, 
|\&„(0)p is independent of n, by (26), so £■„ ~ n^^^. 
Using the connection 



2 2 I ^ " 



M(V0)2 ' 

we can derive a variety of semiclassical sum rules, including 



(28) 



E 



n=narrow 



JO 



TT 



(29) 



where 5 = 2M{Qq) — 2mQ. These are useful in evaluating the heavy-quark 
mass mg, and in calculating the cross section for heavy-quark photoproduction 
using vector-meson-dominance techniques. 

The connection (27) was generalized to higher partial waves by Bell and 
Pasupathy [17], who found 



r=0 



1 

TT 



(2£ + l)!! 



^+V2 d{2t,E^e/h') 



dn 



(30) 



and generalized to include singular potentials by Moxhay & Rosner [18]. 

For a monotonically increasing potential, the semiclassical quantization con- 
dition 



fro , 

dr[2,,{E-V{T))fl^ = {n-\)7: 
connects the shape of the potential to the level density: 



(31) 



dEr, 
dn 



(32) 



Equation (32) is the semiclassical solution to the inverse bound-state problem. 
Although we never applied it to the quarkonium problem, it stimulated us to 
think of ways to reconstruct the interquark potential from the properties of 
the narrow quarkonium levels. 



THE INVERSE BOUND-STATE PROBLEM 



In one space dimension, binding energies and phase shifts uniquely define 
a symmetric potential V{x) = V{—x), for which V{oo) approaches a constant 



E. =3.75GeV 



E.=3.8 GeV 



E. =3.85 GeV 




E =3.9 GeV 









in 

,< - 


T \ 








ol 1 1 ■ ■ 




' 6 






1 — I — 1 1 l_ 








r(GeV 



FIGURE 6. Interquark potentials reconstructed from the masses and leptonic widths of 
V'(3.095) and V'(3.684). The levels of charmonium are indicated on the left-hand side of 
each graph. Those of the T family are shown on the right-hand side of each graph. The 
solid lines denote ^Si levels; dashed lines indicate the levels. (From Thacker, Quigg, 
and Rosner [20].) 



(finite) value. For a "reflectionless" potential (trivial phase shifts), V{x) is an 
algebraic function of the binding energies [19]. For the s-wave inverse problem 
in three dimensions, the central potential is implied by binding energies and 
wave functions at the origin. 

Thacker, Rosner, and I developed a method of successive approximation 
to confining potentials in terms of a sequence of refiectionless potentials that 
support a finite number of bound states [20] . (It is possible to prove interesting 
statements about the convergence of the procedure [21].) Figure 6 shows 
our first attempts to reconstruct potentials from what was known about the 
narrow cc levels, and to use those potentials to predict the properties of bb 
states. In time, we were able to determine potentials separately from the ip 
and T families, and compare them. They agree remarkably well, except at 
the shortest distances, to which the T spectrum has greater sensitivity [22]. 



The inverse-scattering approach is free from assumptions about the short- 
distance and long-distance behavior of the potential. It provided additional 

evidence for flavor independence of the interquark potential, and gave us new 
information on the shape of the potential and where it is determined — for 
0.1 fm^r^l fm. Some important methodological improvements have been 
achieved using techniques of supersymmetric quantum mechanics [23] . 



MESONS WITH BEAUTY AND CHARM 

The next hurrah for quarkonium physics will be the experimental investi- 
gation of the Be family of be bound states. The family is interesting as 
a heavy-heavy system that occupies the region of space between the J/ip and 
the T. Since we know the heavy-quark potential in that region, we should be 
able to calculate the properties of the be states rehably. Unlike the excited cc 
and bb states, all of the be levels below BD threshold are stable against strong 
decays to light flavors (focy^gluons). They cascade by photonic or hadronic 
transitions to the Be ground state. The interest in these states is not merely 
academic. They will soon be discovered and studied at the Tevatron Colhder 
through the decays Be iI^tt, ipiu, . . . 

Estia Eichten and 1 have computed the spectrum of be bound states in a 
number of interquark potentials [24]. We find that the mass of the ground 
state should lie in the range Mb, = 6.258 ±0.020 GeV/c^. The low-lying levels 
in the Buchmiiller-Tye potential [25] are shown in Figure 7. It is noteworthy 
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FIGURE 7. The spectrum of be states in the Buchmiiller-Tye potential (after Eichten & 
Quigg [24]). 
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FIGURE 8. Prominent transitions in the Be spectrum (after Eichten & Quigg [24]). 

that approximately 15 s,p, rf-wave states lie below BD threshold. 

We have also computed the El, Ml, and hadronic transitions between be 
levels. The transitions involving n ~ 2 and n — 1 states are shown in Figure 
8. The narrow widths of these excited states are gathered in Table 2. We find 
that the deep binding of the heavy quarks has a profound influence on decay 
rates of the B^. For example, we estimate that Jb^ ~ 420 MeV '^Sfj^, which 
implies that the purely leptonic decay Be — > rur will be unusually prominent. 

In view of the CDF Collaboration's success in reconstructing the Xc and Xb 
states, I believe that a reasonable — though challenging — experimental goal 
will be to map the eight lowest-lying be states (IS, 2S, 2P) through the 
transitions 2''^Si ^ T'^Si + tttt, 2^80 ^ I^Sq + tttt, Be + 455-MeV 7s, and 
{B* 5e7(72 MeV)) + 353-, 382-, 397-MeV 7s. 

Phenomenological issues raised by the Be family include the systematics of 
spin splittings for the unequal-mass be system and the importance of relativis- 
tic 0{j3'^) corrections. As a third quarkonium system. Be should provide a 

TABLE 2. Total widths of 
low- lying excited states of B^. 



he State Total Width [keV] 



l«Si 


0.135 


2^80 


55 


23Si 


90 


23Po 


79 


2(1+) 


100 


2(1+') 


56 


23P2 


113 




splendid test of a priori calculations from lattice QCD. 

ENVOI 

What have we learned from two decades of quarkonium spectroscopy? The 
first lesson, which underlies all the others, is that nonrelativistic quantum 
mechanics is an appropriate tool for interpreting the quarkonium spectra. 
Using this tool, we have been able to demonstrate by comparing the cc and 
bb systems that the force between quarks is flavor independent, as we expect 
from QCD. Moreover, we have been able to map the interaction between heavy 
quarks in the range 0.1 fm ^ r ^ 1 fm. 

The potential-model approach allows a predictive spectroscopy, including 
calculations of spin splittings. El transition rates, the characteristics of 2S — *• 
IS hadronic transitions, and the properties of wave functions at the origin, 
which are crucial inputs for calculations of quarkonium production in hadronic 
interactions. 

It goes without saying that we have also learned a lot about quantum me- 
chanics! 

What can we hope to learn in the years to come? Among experimental 
goals, we should endeavor to complete the charmonium spectrum: refine our 
knowledge of the r]c and the ^Pj states, confirm the ^Pi (he) level, find the r]'^, 
and search for narrow D-states. It is also desirable to expand our knowledge 
of the T spectrum: locate the rjh and r^^ and the l^Pi (hb) level, and search for 
ID and 2D states. I am optimistic that we shall soon find the and begin 
to explore the be spectrum. 

On the theoretical side, we should be able to refine our understanding of 
relativistic effects and spin splittings. (Information from the spectrum 
should help, in time.) It may be profitable to revisit the coupled-channel effects 
that influence the spectrum near and above the flavor threshold. There may be 
lessons of value for i?-factory experiments here. And finally, we theorists are 
determined to solve QCD and predict the interaction between heavy quarks. 
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